The Droplet State and the Compressibihty Anomaly in Dilute 2D Electron Systems 
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We investigate the space distribution of carrier density and the compressibility of two-dimensional 
(2D) electron systems by using the local density approximation. The strong correlation is simulated 
by the local exchange and correlation energies. A slowly varied disorder potential is applied to 
simulate the disorder effect. We show that the compressibility anomaly observed in 2D systems 
which accompanies the metal-insulator transition can be attributed to the formation of the droplet 
state due to disorder effect at low carrier densities. 



The recent discovery Q of a two-dimensional (2D) 
metal-insulator transition (MIT) has raised the impor- 
tant question concerning the existence of a metallic phase 
in 2D systems. In contrast to the scaling theory of local- 
ization [Q, which predicts that only an insulating phase 
exists in 2D, there is strong experimental evidence |^] 
for metallic-like behavior in many 2D samples. This 
should not be totally surprising because the dominant 
Coulomb interaction in these systems may invalidate the 
non-interacting scaling theory. These intriguing experi- 
ments generate renewed interests in studying the prop- 
erties of low-density 2D electron systems, especially in 
the combined effects of interaction and disorder in such 
systems Most experimental work in the past has 
concentrated on transport measurements. Some recent 
experimental studies |^J^ on thermodynamic properties, 
such as compressibility k in 2D systems, have shed fur- 
ther light on understanding the 2D MIT. It is found 
that the negative 1/k at low densities reaches a minimum 
value at a certain density n, and then increases dramati- 
cally with further decreasing n. Although this surprising 
upturn of 1/k (compressibility anomaly) was observed 
much earlier in a pioneering work by Eisenstein et al. , 
this is the first time that the minimum point in 1/k is 
identified as the critical density for the 2D MIT ||] . On 
the theory side, there are recent efforts [0J^ in addressing 
the interplay between interaction and disorder, and their 
effect in thermodynamic properties. 

In this paper we investigate the space distribution of 
carrier density and the compressibility of 2D electron 
systems by using the local density approximation. The 
strong correlation in such systems is simulated by the lo- 
cal exchange and correlation energies. A slowly varied 
disorder potential is applied to simulate the disorder ef- 
fect. We find that at low average densities electrons form 
a droplet state which is a co-existence phase of high and 
low density regions. We show that the compressibility 
anomaly observed in 2D systems that accompanies the 
metal-insulator transition can be attributed to the for- 
mation of the droplet state . 

To investigate the density distribution of a disordered 
2D electron system, we use the density functional theory. 
The total energy functional reads 

E[n] = Erin] + E,,[n] + Ed[n] + E,[n\ + E,ln]. 

Here Et [n] is the functional of the kinetic energy, E^e [n] 



is the direct Coulomb energy due to the charge inho- 
mogeneity and £'d(n) is the potential energy due to the 
disorder. The strong correlation effect caused by the 
electron-electron interaction is included in the final two 
terms: ^^^[n] is the exchange energy and Ec[n] is the cor- 
relation energy. The ground state density distribution 
can be obtained by minimizing the total energy func- 
tional with respect to the density. 

Using the local density approximation, the total ex- 
change and correlation energies are written as 



E. 



x{c)l 



where e|J(j,)("') exchange (correlation) energy den- 

sity for a homogeneous 2D electron system at a given 
density n, which can be determined by quantum Monte- 
Carlo calculations. In this paper, we use the result from 
Tanatar and Ceperlcy [H, 



- 



1 + aix 



1 + aix + a2X^ + flacc^ 



where x = l/(7rn)^/*. The energy unit is IRy* 
m* /2e*'^Ti'^ ^ and the parameters oq — —0.3568, ai 
1.13, a2 = 0.9052, 03 = 0.4165. 

The kinetic energy functional can be written as 

EtH = J dx^^l(x) (-V2) V',(x), 



where the sum is over all occupied quasi-particle energy 
levels, and n(x) — IV'j(x)P- To further simplify the 
calculation, we make an approximation to the kinetic en- 
ergy so that it can be written in the form of a density 
bmctional [0, 



Erin] 



(he 



7rn(x) 



1 |Vn(x)p 
4 n(x) 



The first term provides the local density approximation 
for the kinetic energy, while the second term includes 
the effect of the density gradient. The approximation 
provides enough accuracy for this class of calculations. 
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The energy functional for the disorder potential Vd(x) 
can be written as 



EM 



(ixT/d(x)n(x) . 



In a real system, a disorder potential may be slowly vary- 
ing and has correlation between different positions. To 
simulate the situation, we assume the correlation for the 
disorder follows the simple behavior. 



(Fd(x)^^d(x')) =l^,'exp 



where Vg is the amplitude of the potential fluctuation, 
and ^ is the correlation length of the disorder. ^ is 
roughly the average size of valleys in a disorder land- 
scape. 

In summary, the total energy functional is of the form 



E[n\ = / dx 



^n(x)^ + 7 



1 |Vn(x)|^ 



dx 



,n(x)n(x') 



4 n(x) 

+V^d(x)n(x) -f eS(n)n(x) + e°(n)n(x)] . 

Thus, the local density approximation converts the 
strong-interacting problem to a single particle problem 
with a self-consistently determined potential. The den- 
sity distribution of the ground state can be obtained by 
minimizing the energy functional under the constraint 
of a constant total electron number. We introduce the 
variable x with n(x) = A^x(x)^/ / '^x'x(x')^, where N 
is the total number of the electrons in the system. The 
constraint for the constant total electron number is auto- 
matically satisfied with the new variable. We can get the 
minimized energy functional by using steepest descent 
method with iterations |llt]. 



X 



'"+i(x)=x"(x)-7 



lx=X'"(x) ' 



where 7 is the iteration constant which is chosen so that 
the interaction is convergent, and 



2N 



Sx J dxx(x)^ 



[-V' + V,ff[n]-Eo] x(x). 



where 



Veffin] = 2nn -f Vd{x) + — [eS[n]n + eO[n]n] 



and 



En — 



/dxx(x) [-V'' + Veff[n]] x(x) 
/ dxx(x)2 



The calculation is carried out in a 128 x 128 discrete 
space. The size of the system is set as L = 256o|j, where 
a*g is the effective Bohr's radius, a|j = efi' jme^^ with e 
being the dielectric constant and m the effective mass of 



an electron. The periodic boundary condition and the 
Ewald sum for the Coulomb interaction are applied to 
minimize the finite size effect. The electron density is 
adjusted by changing the total electron number iV. The 
density distribution and the total energy of the system 
are calculated for different densities, and the chemical 
potential is calculated by using the formula, 

^(iV) = E{N + \) - E{N) . 

The compressibility of the system can be calculated by 

1 _ iV^ 9/1 

where S is the total area of the system. 

Figure |^ shows the density distribution of the system. 
It can be clearly seen that the electrons form some high 
density regions, while the density of other regions are 
essentially zero. Depending on the average density of 
the system, the high density regions may connect each 
other together (r^ — 10), or form some isolated regions 
[vg = 18). There exists a certain density (r^ = 14) where 
the high density regions starts to percolate through the 
system, and form a conducting channel. The calcula- 
tion clearly demonstrates the idea of our earlier theory 
i.e., the metal-insulator transition observed in the 
2D electron systems is the percolation transition of the 
electron density. 

The electron-electron interaction is important for the 
conducting behavior of a dilute electron system in the 
sense that it makes the density distribution more ex- 
tended because of the Coulomb repulsion. Figure]^ shows 
the density distribution for the free electron gas with the 
same density as in Fig. 1(b) by turning off the electron- 
electron interaction. The system only forms some iso- 
lated high density regions at the disorder valleys, while 
the density distribution of the corresponding interacting 
system (Fig. 1(b)) is quite extensive at the same density. 
At a given disorder strength, the critical density for the 
free electron gas is much higher than its interacting coun- 
terpart. 

Figure || shows the compressibility of the systems. To 
compare with the experiments 0,^, we calculate 6n/6N, 
which is the direct measured quantity in the experiments. 
It is well known that the compressibility of a uniform 
electron gas is negative in the low density region due to 
the effect of the exchange and correlation energies, as 
shown by the solid line in Fig. 3(a). However, when the 
disorder is present, the behavior changes greatly. In the 
low density, the electrons tend to occupy the valleys of 
the disorder landscape, and the local density, instead of 
the average density, determines the compressibility of the 
system. On the other hand, at higher densities, all of the 
valleys are filled, one can expect the compressibility of 
the system to resume the behavior of a uniform electron 
gas. We have a non-monotonic behavior for 6^/6N, as 
shown by the dots in Fig. 3(a), which are in good agree- 
ment with the experimental measurement [^j6|. Com- 
paring with FigJ|, we find that the turning point of the 
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FIG. 2. The density distribution for the free electron gas on 
the same disorder landscape as the Fig.^ at density = 14. 
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FIG. 1. The density distributions for the different electron 
densities. We use the contour plot for the local density pa- 
rameter Vs = 1/ ^nn. The density of the white area decreases 
rapidly to zero. We take Vs = 0.2Ry* , ^ = 0.2L. 



FIG. 3. Sn/SN as a function of the electron density. Solid 
lines are for the uniform electron gas, squares are the data 
points for the disordered system, and the dashed-lines are 
the results from the theory discussed in the text. Two kinds 
of the disorder are used in the calculation: (a) The same 
disorder landscape as in Figure |l[ (b) Off-plane charge im- 



purities with d = Wag, rn = 2.5 x 10 /a 
eters in the dashed-lines: (a) no = 10" 



The param- 
a = 1.5; (b) 



no — 0.5 X 10 /a*g, a — 2.3. is the total number of 



the electrons in the simulation box. 
Tb = 18; = 50 corresponds to r. 
sponds to = 10. 



N 



s to 

= 30 corresponds to 
14; N = 100 corre- 
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compressibility {N « 50, Ts « 14) coincides with the 
percolation threshold of the system. At low densities, 
the data points in the plot show strong fluctuation, indi- 
cating the effect of the local fluctuation of the disorder 
potential. 

The behavior can be understood by a simple theory. 
Following the definition, the chemical potential fi is the 
energy needed to add an electron into the system, 

fi{N) ^ E{N+l) - E{N) 
«eo [neff{N+l)] {N- 
_ _S_ 
~ SN 

where we suppose that the electron energy is determined 
by the local density of the electrons. E{N) is the to- 
tal energy of the system, eo{n) is the energy per elec- 
tron for the uniform electron gas, and neff is effective 
local density. For the inhomogeneous system as shown 
in Fig.|l|, the effective local density can be estimated by 
neff{n) ~ n/ fin), where /(n) is the fraction of the high 
density region. After some algebra, we have 



l)-eo [neff{N)]N 



{eo Wff{N)]N}. 



5fJL 



1 - 



din/ 
dlnn 



- £o 



din/ 
dlnrt' 



1 - 



din/ 
dlnn 



where /ig is the chemical potential for a uniform electron 
gas. In the low density limit, f{n) — > 0, the local effec- 
tive density is greatly different from the average density 
of the system. As a consequence, the density dependence 
of the chemical potential, 5^/5n, changes greatly. In gen- 
eral, supposing /(n) n" in the low density limit, the 
analysis shows that Sfi/Sn will have a non-monotonic be- 
havior if a > 1. The behavior of f{n) is determined by 
the local disorder potential profile. In a 2D system, the 
infinite harmonic potential has f{n) n. So the require- 
ment a > 1 is equivalent to the condition that the local 
disorder potential has a weaker confinement effect than 
the harmonic potential. The condition can be easily sat- 
isfied in a typical experimental system. For instance, the 
coulomb potential has f{n) ^ n^. In Fig.||, we use the 
above equation for Sfj,/Sn with the following relation of 
/(n) to fit the data, 



1 



1 + ^ 



This form has a correct behavior in the high density limit, 
/ — > 1, and the low density behavior is controlled by a. 
By carefully choosing the values for a and no, we obtain 
a good agreement with our numerical data as shown in 
the dashed line in Fig. 3(a). 

To demonstrate the effect of a different disorder poten- 
tial, we use off-plane charge impurity potential in calcu- 
lating Figure ||(b), i.e., 

Vd(x) = ^ ^ 



E 



where d is the distance between the electron and the 
impurity planes, and the impurities are randomly dis- 
tributed with a density n^. This potential gives a similar 
density profile as shown in Figure 1. As expected, this 
form of potential has a larger value of a as given in the 
figure caption. 

To conclude, we have studied the electron space dis- 
tribution and the compressibility of disordered dilute 2D 
electron systems by using the local density approxima- 
tion. Electron distribution confirms the formation of the 
droplet state that consists of high and low density re- 
gions. Our calculated compressibility is in good agree- 
ment with the experimentally observed behavior showing 
unexpected anomaly at low densities. The turning point 
of the compressibility happens around the percolation 
threshold. Our theory based on the droplet state provides 
a good understanding of the compressibility anomaly. 

This work is supported by DOE. 



[2] 
[3] 

[4] 
[5] 
[6] 

[7] 
[8] 

[9] 

[10] 

[11] 



S. V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, 
V. M. Pudalov, and M. D'lorio, Phys. Rev. B50, 8039 
(1994); S. V. Kravchenko, W. E. Mason, G. E. Bowker, 
J. E. Furneaux, V. M. Pudalov, and M. D'lorio, Phys. 
Rev. B51, 7038 (1995); S. V. Kravchenko, D. Simonian, 
M. P. Sarachik, W. Mason, and J. E. Furneaux, Phys. 
Rev. Lett. 77, 4938 (1996). 

E. Abrahams, P.W. Anderson, D.C. Licciardello, and 
T.V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
For a review, see E. Abrahams, S. V. Kravchenko, 
and M.P. Sarachik, :ond- mat/0006054 and references 
therein. 

C. Dultz and H.W. Jiang, Phys. Rev. Lett. 84, 4689 
(2000). 

S. Hani, A. Yacoby, D. Mahalu, and H. Shtrikman, Phys. 
Rev. Lett. 84, 3133 (2000). 

J.P. Eisenstein, L.N. Pfeiffer, and K.W. West, Pliys. Rev. 
Lett. 68, 674 (1992); J.P. Eisenstein, L.N. Pfeiffer, and 
K.W. West, Phys. Rev. B50, 1760 (1994). 
Q. Si and CM. Varma, Phys. Rev. Lett. 81, 4951 (1998) 

A. A. Pastor and V. Dobrosavljevic, Phys. Rev. Lett. 83 
4642 (1999). 

S. He and X.C. Xie, Phys. Rev. Lett. 80, 3324 (1998); J 
Shi, S. He and X.C. Xie, Phys. Rev. B60, R13950 (1999) 

B. Tanatar and D.M. Ceperley, Phys. Rev. B39, 5005 
(1989); S.T. Chui and B. Tanatar, Phys. Rev. Lett. 74 
458 (1995). 

Y. Wang, J. Wang, H. Guo and E. Zarembar, Phys. Rev 
B52, 2738 (1995), and references therein. 



<P 



4 



